Critical behaviour of the Id annihilation fission process 2A — > 0, 2A — > 
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Numerical simulations and cluster mean-field approximations with coherent anomaly extrapolation 
show that the critical line of the Id annihilation fission process is separated into two regions. In both 
the small and high diffusion cases the critical behavior is different from the well known universality 
classes of non-equilibrium phase transitions to absorbing states. The high diffusion region seems 
to be well described by the cyclically coupled directed percolation and annihilating random walk. 
Spreading exponents show non-universal behavior. 



PACS numbers: 05.70.Ln, 82.20.Wt. 

Non-equilibrium phase transition may take place even 
in one dimensional systems. However the ordered phase 
lacks fluctuations that could destroy the state. If the sys- 
tem has evolved into that state it will be trapped there. 
We call this transition to an absorbing state. First order 
transition among these is very rare, there have been only 
a few system found that exhibit a discontinuous jump 
from the active to the absorbing phase Jl]-|5| . 

Continuous phase transitions have been found to be- 
long to a few universality classes, the most robust of them 
is the directed percolation (DP) class. According to the 
hypothesis of ||J^] all continuous phase transitions to a 
single absorbing state in homogeneous systems with short 
ranged interactions belong to this class provided there 
is no additional symmetry and quenched randomness 
present. The most prominent system exhibiting phase 
transition of this class is the branching and annihilating 
random walk with one offspring (BARW). Furthermore 
systems with infinitely many absorbing states were also 
found to belong to this class 

An important exception from the DP class is the so- 
called parity-conserving (PC) class where particles fol- 
low a branching and annihilating random walk with two 
offsprings (BARW2) A -> 3A, 2A -> that conserves 
the parity of their number ]To|~[r^|. Particles following 
BARW2 may also appear as kinks between ordered do- 
mains in systems exhibiting two absorbing states Q [l8| . 
However it was realized that the BARW2 dynamics alone 
is not sufficient condition for a transition in this class but 
an exact Z2 symmetry between the two absorbing states 
is necessary, too p?| , pj| -p2| . 

Recently a study on the annihilation fission (AF) pro- 
cess 2 A -> 0, 2 A -> 3 A H| suggested that neither the 
BARW2 dynamics nor the Z% symmetry but simply the 
occurrence of two absorbing states can result in a PC 
class transition. This model was introduced first by J2J] 
in the context of a renormalization group analysis of the 
corresponding bosonic field theory. The main idea was 
to interpolate between systems with real and imaginary 
noise. These studies predicted a non-DP class transition, 
but they could not tell to which universality class this 
transition really belongs. 



A fermionic version of the AF process introduced in 
p3| is controlled by two parameters, namely the proba- 
bility of pair annihilation p and the probability of particle 
diffusion d. The dynamical rules are 



AA0, 0AA -> AAA 
AA^00 
A0 0A 



with rate(l - p){l - d)/2 
with ratep(l — d) (1) 
with rated . 



Carlon et al. suggested that the phase diagram can be 
separated into two regions. For low values of d (less than 
approximately 0.3), they found a continuous phase tran- 
sition belonging to the PC universality class. For large 
values of d, however, a first order transition was reported. 
This claim is based on mean-field, pair mean-field and 
density matrix renormalization group method (DMRG) 
|p5| calculations. An even more recent preliminary nu- 
merical study H^] found non-PC class like critical expo- 
nents and posed the question if there is a new type of 
non-equilibrium phase transition occurring here. 

In this paper I report mored detailed numerical results: 
Simulations and generalized mean-field (GMF) approxi- 
mations with coherent anomaly (CAM) extrapolations 
that give numerical evidence of a rich phase diagram with 
two different kinds of new universality classes in the low 
d and high d regions. 

For simulations a parallel update version of the AF pro- 
cess was used since the model can effectively be mapped 
on a massively parallel processor ring |27j]. To avoid 
collisions by simultaneous updates the annihilation and 
exchange steps are done on two sub-lattices, while cre- 
ation is performed on three sub-lattices (with a rate 
double as in eqjj]). The critical point has been deter- 
mined for d = 0.05,0.1,0.2,0.5,0.9 by following the de- 
cay of the particles from a uniform, random initial state 
(p(t) oc t~ 5 ). The local slopes curve of the density 



Seff(t) = 



-\n[p(t)/p(t/m)} 
ln(m) 



(2) 



(where we use m — 8 usually) at the critical point goes to 
exponent S by a straight line, while in sub(super)-critical 
cases it veers down(up) respectively. As Fig. |l| shows at 
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d = 0.2 and p = 0.24802 this quantity scales without any 
relevant corrections (i.e. the straight line is horizontal) 
and in the 1/t ->• limit goes to 5 = 0.268(2). Very 




0.0001 0.0002 0.0003 0.0004 

1/t 



FIG. 1. Local slopes of the particle density decay in the 
parallel update model at d — 0.2. Different curves correspond 
top = 0.2481, 0.24807, 0.24805, 0.24802, 0.248, 0.24795 (from 
bottom to top). The simulations were performed on a ring of 
size L = 24000 up to 4 x 10 5 MCS. Throughout the whole 
paper t is measured in units of Monte-Carlo sweeps. 

similar results has been obtained for d — 0.05 and d = 0.1 
showing that the exponent 8 is close to the corresponding 
value of the PC class (0.285(10) Q), but appears to be 
significantly smaller. This 8 is in agreement within error 
margin with the value (<~ 0.272(18)) one can obtain from 
DMGR calculations |2^] using scaling laws S = &/v\\ = 
0/"±/z. 

For stronger diffusion rates d = 0.5 and d = 0.9 
one can observe slower density decays with exponent 
8 = 0.215(10) (Fig. |). Note that I checked, that chang- 
ing the abscissa on the local slopes graphs from 1/t to 
hx(t)/t does not eliminate the differences between small 
and large d behaviors. Also by assuming strong correc- 
tion to scaling I could not get coherent estimates for (3 
and 8 that would be in agreement with DP or PC classes. 
This value of 8 is in agreement with that of the cou- 
pled DP+annihilation random walk suggested by [ESj. 
This model was proposed to explain the space-time be- 
havior of the AF with the help of a multi-component 
model, where A particles (corresponding to pairs in AF) 
perform BARW process: A —* 2A, A — > and occa- 
sionally create B particles (corresponding to lonely par- 
ticles in AF) who follow random walk + annihilation : 
2B — > A (ARW). The space-time evolution picture of 
this model looks similar to that of AF (and very differ- 
ent from that of BARW2) but the numerical simulations 
predicted somewhat different exponents as those of the 
AF p3| , p6| , p9|j3^ ] . By approaching d — 1 corrections to 
scaling are getting stronger and this might leaded to the 
conclusion of pj|| that for d > 0.3 the transition is of first 
order. The present study can not confirm those predic- 
tions based on extrapolations to DMRG results of small 



system sizes (L < 100). The first order transition for 
d > 0.3 appears to be unlikely also because the mean-field 
approximation of the model exhibits continuous phase 
transition (with (3mf = 1) and that corresponds to the 
strong diffusion limit without fluctuations. 
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FIG. 2. The same as Fig. [l] at d = 0.5. Different curves 
correspond to p = 0.2796, 0.27958, 0.27955, 0.2795 (from bot- 
tom to top). 

The differences from PC class and from first order tran- 
sition can be seen even more clearly by direct measure- 
ments of the order parameter (3 exponent (p(oo) oc e' 3 ). 
By looking at the effective exponent defined as 

o ( \ lnp(e») -lnp(ej-i) 

Ine.-lne^ ' (3) 

on Fig. H one can observe two regions again: for d < 0.3 
the f3 e ff tends to 0.58(2) as e = (p c — p) — ► 0, while for 
d > 0.3 it goes to ~ 0.4(2). 

For small diffusion rates the value 0.58 is very different 
from that of the PC class (0.94(10)) @||) and from 
that of DP class (0.2765(1)) @, but is in agreement with 
preliminary simulation results of [ ^6| and po|] . For strong 
diffusion it coincides fairly well with the value 0.38(6) 
that was found for the coupled DP+ARW model |2S[] . 

The results for the order parameter exponent has been 
verified by generalized mean-field (GMF) approximations 
||3"i"| , j32"f with coherent anomaly extrapolation (CAM) [j33[ . 
This method has been proven to give precise estimates 
for the DP @j and PC ^classes. The details of the 
method are described in ]22fl and a forthcoming paper 
p5[ will discuss the results and show estimates for other 
exponents as well. 

For technical reasons a variant of the AF model was 
investigated where the creation is symmetric: A0A — » 
AAA. This does not change the essence since the 
AA0/0AA — * AAA process can be decomposed into a 
A0 — > 0A exchange and a symmetric creation. One can 
easily check that the mean field approximation and result 
of it is the same 

p(oo) = (l-3p)/(l-p) (4) 
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FIG. 3. Order parameter exponent results of the AF model 
for d = 0.9 (circles), d = 0.5 (squares), d = 0.2 (diamonds), 
d = 0.1 (triangles), d = 0.05 (stars). The simulations were 
performed on a ring of L — 24000 sites such that averaging 
was done following the steady state is reached for 1000 sam- 
ples 

Higher order approximations for N — 2, 3, ..7 converge 
to the simulation results of this variant (Fig. |J) and the 
phase diagram qualitatively agrees with that of ref . |23] . 
The leading order singularity at critical point remains 
mean- field type i.e. (3gmf — 1. Furthermore dynami- 
cal MC simulation of this variant gives the same critical 
decay behavior as the original AF model. According to 
CAM the amplitudes of these singularities a(N) scale in 
such a way that 

a{N)o,{ Pc {N)-p c f-^T (5) 

the exponent of true singular behavior can be estimated. 
The CAM extrapolation has been applied for the high- 
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FIG. 4. GMF approximations of the particle density in the 
symmetrized AF model at d = 0.5 for N = 1,2, ...7 (from 
right to left curves). 

est levels of approximations [N — 4, 5, 6, 7) and resulted 



TABLE I. Summary of results. The non-universal critical 
parameter p c of the parallel model is shown here. 



d 


0.05 


0.1 


0.2 


0.5 


0.9 


Pc 


0.25078 


0.24889 


0.24802 


0.27955 


0.4324 


PCAM 




0.58(6) 


0.58(2) 


0.42(4) 






0.57(2) 


0.58(1) 


0.58(1) 


0.4(2) 


0.39(2) 


S 


0.273(2) 


0.275(4) 


0.268(2) 


0.21(1) 


0.20(1) 


V 


0.10(2) 




0.14(1) 


0.23(2) 


0.48(1) 


51 


0.004(6) 




0.004(6) 


0.008(9) 


0.01(1) 



in [3cam values is agreement with those of simulations 
(see summarizing Table |j). Note, that the CAM extrap- 
olation does not work for d > 0.5 and d < 0.1, because 
corrections to scaling are getting strong, and it is very 
difficult to find the steady state solutions of the GMF 
equations for N > 7. 

Time dependent simulations from a single active seed 
for these models have been proven to be a very efficient 
method |3i| since the slowly spreading clusters do not 
exceed the allowed system sizes hence do not feel finite 
size effects. Very precise estimates for the cluster "mass" 
exponent N(t) cx f and the cluster survival probabil- 
ity exponent P(t) oc t~ 5 ' have been obtained in DP and 
PC classes (see references in |37|). However example in 
the case of systems with infinitely many absorbing states 
these exponents seem to depend on initial conditions 
pLpl. A rigorous explanation is still missing. The simu- 
lation results for the cluster mass (Fig. |[|) show strong 
dependence of i] on d. For d = 0.5 and d = 0.9, where 
(3 and 6 are roughly the same it changes from 0.23(2) to 
0.48(1). The survival probability exponent 8i is roughly 
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FIG. 5. Local slopes of the particle number in seed sim- 
ulations in the parallel update model at d = 0.05, d — 0.2, 
d — 0.5 and d — 0.9 (from bottom to top). The simulations 
were performed on a lattice of size L — 24000 up to 10 5 MCS. 

zero, expressing the fact that a single particle survives 
with finite probability (see table ^) . 

The continuously changing exponent r\ is in agreement 
with the preliminary observations of pq ] , however (3 and 
the 6 seem to be constant in the region: 0.05 < d < 0.2 
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and 0.5 < d < 0.9, suggesting two different universality 
classes. 

The non-universal behavior of the cluster exponent 
may not be so surprising since it was also observed in 
the d = limit of the AF in the pair contact process 

Si' 

In conclusion numerical simulations and analytical 
GMF+CAM calculations suggest two new universality 
classes in the AF model along the phase transition line. 
For small diffusion single particles cannot escape clus- 
ters where pair annihilation and creation dominates. 
For larger diffusions these single particles will wander 
among pair creation-annihilation clusters and the cou- 
pled BARW + ARW description of §§ will be valid. 
These new universal behaviors may sound uneasy follow- 
ing two decades when the DP hypothesis has been found 
to describe almost every continuous transition to absorb- 
ing state, but do not contradict to it at all. In the AF 
model the absorbing state is not singlet, there are two 
(non symmetric) absorbing states where the system can 
evolve with equal probabilities. One of them is not com- 
pletely frozen, but a single particle diffuses in it. My re- 
sults do not support the conclusions of Q for a PC class 
transition without BARW2 process and Z2 symmetry nor 
the first order transition for large d, but the exponents 
do not contradict to those of p3|] within numerical preci- 
sion. It is likely that the two exponent ratios determined 
by |2J| arc close to those of PC class accidentally. The 
cluster exponent 77 behaves non-universally like in the 
pair contact process, the d = limit of the AF. Deeper 
understanding of this model that would reveal hidden 
symmetries responsible for the strange critical behavior 
of AF would be highly desirable. 
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